Energy levels of a two-dimensional hydrogen atom with spin-orbit Rashba interaction 
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Electronic bound states around charged impurities in two-dimensional systems with structural 
inversion asymmetry can be described in terms of a two-dimensional hydrogen atom in the presence 
of a Rashba spin-orbit interaction. Here, the energy levels of the bound electron are evaluated 
numerically as a function of the spin-orbit interaction, and analytic expressions for the weak and 
strong spin-orbit coupling limits are compared with the numerical results. It is found that, besides 
the level splitting due to the lack of inversion symmetry, the energy levels are lowered for sufficiently 
strong spin-orbit coupling, indicating that the electron gets more tightly bound to the ion as the 
spin-orbit interaction increases. Similarities and differences with respect to the two-dimensional 
Frohlich polaron with Rashba coupling are discussed. 
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The two-dimensional (2D) hydrogen atom, i. e., an 
electron constrained to move in a plane and subjected 
to an attractive Coulomb potential j L S i j m i is a theoreti- 
cal construction which, besides being of interest in itself, 
has also important physical realizations. It can describe 
indeed the effect of a charged impurity in 2D systems 
such as quantum wells and surface states, or in extremely 
anisotropic three-dimensional crystals, 1 as well as exci- 
tons in semiconductor 2D heterostructures£ 

The spin-orbit (SO) interaction, arising from the struc- 
tural and/or bulk inversion asymmetries, characterizes 
several of the above mentioned low-dimensional systems^ 
and gives rise to energy level splittings ranging from a few 
to hundreds of meV, depending on the material char- 
acteristics (see for example Ref.Q)- Furthermore, the 
possibility of tuning the SO interaction in semiconduc- 
tor quantum wells by means of external applied voltages 
represents the key feature for application in spintronics. 
Given this situation, it becomes therefore natural to as- 
sess how the properties of a 2D hydrogen atom are af- 
fected by the SO interaction. 

Several studies have already been devoted to the ef- 
fect of the SO coupling in electrons interacting with 
central potentials, such as those describing hard-wall or 
parabolic quantum dots i 8 ' 9 ' 10 ! 11 ! 12 However, despite its 
potential interest for SO coupled low-dimensional sys- 
tems, the specific 2D Coulomb problem appears to have 
been only marginally considered in the literature^ In 
this Brief Report, the 2D Coulomb problem is numer- 
ically solved for an electron interacting with a Rashba 
potential, that is the SO coupling arising from structural 
inversion asymmetry in the direction perpendicular to 
the 2D planed It is found that the Rashba interaction 
removes partially the initial degeneracy of the 2D hydro- 
gen atom, and the resulting energy levels are two-fold de- 
generate due to the time-reversal invariance of the model. 
Furthermore, it is shown that the SO interaction renders 
the electron more tightly bound to the ion, confirming 
a general trend observed for other central potentials and 
for 2D electrons coupled to phonons. 

The Hamiltonian for a 2D hydrogen atom in the pres- 



ence of a Rashba SO potential is as follows (h = 1) 
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where p q = —id/ dq is the electron momentum operator 
(q = x, y), p 2 = p x + p y 2 , m e is the electron mass, and 
a x and a y are Pauli matrices. The last term of Eq.([T]) de- 
scribes the Rashba SO interaction with coupling param- 
eter 7. For 7^0 but zero Coulomb interaction [Z = 0), 
Eq.(fT]) is easily diagonalized in the momentum space, and 
the resulting energy dispersion of the free electron is com- 
posed of two branches Ek,± = (k ± kn) 2 /2m e — Er, 
where kfj = m e 7 is the Rashba momentum. In the 
ground state, the electron has energy Ek R - = —Er, 
where Er = k R /2m e = m e j 2 /2. 

In the presence of the Coulomb interaction, it is con- 
venient to rewrite Eq.JTJ) in polar coordinates: 
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is the free electron Hamiltonian. Equation ^ commutes 
with the z-projection of the total angular momentum 
J z = L z + <j z /2, where L z = —id/d<f), so that the eigen- 
functions of ([2]) can be chosen to be simultaneously eigen- 
functions of J z . Since H in polar coordinates allows for 
separation of variables, its eigenfunctions have therefore 
the following form£&±&±i^ 



i(j-l/2)0 
t (3+1/2)* 



(4) 



where j = ±1/2, ±3/2, . . . are the eigenvalues of J z . The 
lack of spatial inversion symmetry induced by the pres- 
ence of the Rashba interaction lowers the symmetry of 



2 



H when 7^0. As shown below, this will induce a 
splitting of the energy levels compared to the case when 
7 = 0. Note however that H commutes with the time- 
reversal operator K = ia y C, where C is the operation of 
complex conjugation, so that and its Kramer conju- 
gate Ktyj have the same energy. This implies that, since 
J z ^fj = j^j and J z Ktyj = —jK^j, the energy spectrum 
of H is invariant under the change j — » —j. 

For bound states, the Schrodinger equation H^Sfj = 
Ety j is rewritten by introducing q$ = —2m e E and the 



dimcnsionlcss radial variable p 
one therefore finds 



2q r. By using Eq.@ 
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2m e Ze 2 . The first terms in Eqs.© and © 
represent the differential equations for the radial wave 
function of the usual 2D Coulomb problem (i. e., with- 
out SO coupling) with quantum numbers j — 1/2 and 
j + 1/2, respectively!^ Apart from a normalization con- 
stant, their solutions are of the form Rnj±i/2(p) = 
cxp(-p/2)p^ ±1 ^L 2 I !/ ±1/21 (p) where N = 0,1,2,... is 

the radial quantum number and L 2 ^^ 1 ^ 2 ^ (p) arc Laguerre 
polynomials The corresponding energy levels are 

r?/4 



(N+ | j ± 
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where r] = 2m e Z 2 e 4 = A 2 /2m e . By introducing the prin- 
cipal quantum number n = N+ | j ± 1/2 |= 0, 1, 2 . . ., 
with | j ± 1/2 |< n, one infers that each level with fixed 
n has energy —r\j (2n + l) 2 with degeneracy 2(2n + 1). 

For nonzero SO coupling, it is natural to expand the 
radial functions /• in terms of Rnj±i/2(p)- By keeping 
in mind that the total wave function ^ j must be also 
cigenfunction of J z , one has: 

00 

ff{p)=e-^p\^Y. A %A j±m (P)- (8) 

By substituting the above expansions in Eqs.© and 
and by making use of the properties of the Laguerre 
polynomials^ one arrives at the following iterative sys- 
tem of equations 
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FIG. 1: Energy levels obtained from a numerical solution of 
Eqs.Q and (|10|l (solid lines) as a function of the Rashba en- 
ergy Er. The dashed lines are the analytic results of Eq. (|13p 
for the weak SO limit. All energy levels are shifted with re- 
spect to the ground state of the free electron with SO interac- 
tion [Eg , <|14p ] . The different levels are labeled by the principal 
quantum number n, reported in the left axis, and by the ra- 
dial quantum number N and the total angular momentum in 
the z direction j reported in parentheses. 
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The values of qo, and so the energy levels E = — (7 2 /2m e , 
satisfying Eqs.© and (fT0|) can be easily obtained an- 
alytically in the limit of weak SO coupling. It suffices 
to recognize that decoupling Eqs.© and (p~0|) leads to 
two iterative equations of the form a± A^ . + b± A^_ 2 ■ + 
c±A^ +2 j = 0, whose solutions in the weak SO limit arc 
determined simply by the condition a± = 0, since b± and 
c± are both of order (fcfj/A) 2 . Up to order (&r/A) 2 the 
coefficients a± are given by 
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so that the energy levels £wj±i/2 of the weak SO inter- 
acting case are 



E 



N,j±l/2 



E N,j±i/2 ± 2jER - E R , 
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where E° N j± 1 / 2 ^ s the energy spectrum for zero SO inter- 
action given in Eq.((7|). From Eq. (fT3"|) one sees therefore 
that the 2(2n + l)-fold degeneracy for zero SO coupling 
is lifted when 7^0 and that each level is splittcd into 
2n + 1 levels, each two-fold degenerate. The remaining 
degeneracy is due to the time-reversal invariancc of H 
and can be removed by adding a time-reversal breaking 
term in the Hamiltonian such as a magnetic field. Note 
also that Eq. (fT3|) could be obtained by making use of the 
method described in Ref . [l~5T| . 

A comparison between the weak SO coupling result 
(fTB"]) (dashed lines) and the energy levels obtained by a 
numerical solution of Eqs.© and (|10j) (solid lines) is re- 
ported in Figfl] for a few energy levels. The principal 
quantum number values n are reported in the left axis, 
while the radial and total angular momentum quantum 
numbers TV and j are indicated in parentheses. What is 
plotted in Figfl] is actually the quantity 



£iV,j±l/2 — ^N,j±l/2 + Er, 



(14) 



that is the energy spectrum shifted with respect to the 
ground state — En of the free electron coupled to the 
SO potential. As it is shown in the figure, the ground 
state (identified by quantum numbers n = 0, TV = 0, 
and j = ±1/2) has its energy lowered by the SO interac- 
tion, demonstrating that the electron gets more tightly 
bounded as En increase. This holds true also for the 
higher energy levels which, besides being splitted by the 
Rashba interaction, have their energies lowered for suf- 
ficiently large En values, as it apparent for most of the 
levels plotted in Figfl] For states like n = 1 (0, ±1/2) 
and n = 2 (2, ±1/2) one needs En/rj > 2 before reaching 
net energy levels lower than the zero SO limit. This is 
of course unattainable since, for an unscreened charged 
impurity, r\ is of the order of one Ry, while the maximum 
value of En to date is of about 0.2 eV J£ The relevant val- 
ues of E r /r] are therefore lower than about 0.01 for which, 
however, the perturbative result (|14[) is in quantitative 
good agreement with the numerical solution plotted in 
Figfl] 

Let us now discuss the relevance of the above results 
with respect to a different but somewhat related problem: 
that of a 2D electron strongly coupled to longitudinal 
optical phonons in the presence of a SO Rashba poten- 
tial (Frohlich-Rashba model) ill To this end it is useful 
to compare the exact numerical solutions of Eqs.® and 
(fi"0)) with a variational calculation of the ground state 
energy Eq obtained from the following electron trial wave 
function: 



Mbr) 
Ji(br) 



(15) 



where a and b are variational parameters, Jo and J\ are 
Bessel functions, and A is a normalization constant. The 
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FIG. 2: Comparison between the exact ground state energy 
obtained from a numerical solution of Eqs.Q and (|10p (solid 
line) and the variational calculation with the ansatz wave 
function (|15[) (dashed line). Inset, the same results plotted 
as a function of ln 2 (ER/r]). 



above form for 4 , o( r I </>) was introduced in Ref. [17| in or- 
der to find an upper bound for the ground state of a 
the Frohlich-Rashba model. By using Eq. lfil)]) , the en- 
ergy functional T = (^0 I {H + En) \ ^0) can be found 
analytically 
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(16) 

where K and E are complete elliptic integrals of the first 
and second kind, respectively. Numerical minimization of 
Eq. (|16j) with respect to a and b provides an upper bound 
£ for the ground state energy. This is plotted in Figf2] 
by the dashed line and compared with the exact ground 
state energy (solid line) obtained from Eqs.© and (|10[) . 
It is seen that the simple ansatz (|15[) reproduces fairly 
well the lowering of the ground state energy as En in- 
creases. In particular, by expanding Eq. (|16[) for small 
values of En compared to ij = A 2 /2m e it turns out that 
£0 ~ — rj — En, which coincides with Eqs. lfTB"]) and (fT4|) 
for n = 0, N = and j = 1/2. For very large values of 
En/rj, Eq. (fl6]) has the limiting form 
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where e is the Neper number and should not be con- 
fused with the electron charge. Equation (jTTJ) is min- 
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imizcd by setting b — kn and, within logarithmic ac- 
curacy, a ~ A ln(4fc#/Ae 2 ). Hence, the corresponding 
asymptotic upper bound for the ground state energy re- 
duces to: 
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indicating that for Er —> oo the ground state energy gets 
indefinitely lowered by following a squared logarithmic 
dependence on Er. This result is confirmed in the inset 
of Fig|5J where both the exact result (solid line) and the 
numerical minimization of (|16p (dashed line) reduce to 
straight lines when plotted as a function of In 2 (£72/77). 

The functional dependence of the ground state energy 
on the SO coupling shown in Eg. (TX5|) was originally ob- 
tained by a different method in Ref.[l 2j where, however, 
a cutoff parameter was introduced to prevent a diverging 
result. As further noted in Ref. [T3| , a squared logarithmic 
behavior characterizes also the ground state energy of the 



three-dimensional (3D) hydrogen atom in an extremely 
strong magnetic field iJ , 18 ' 19 supporting the interpreta- 
tion that a 2D electron in the presence of a strong Rashba 
SO interaction behaves effectively as a one-dimensional 
(ID) particle. The correspondence between ID-like be- 
havior and strong Rashba interaction has been notices 
also for bound states of 2D electrons in short range cen- 
tral potentials ; 12 ' 20 as well as for 2D electrons coupled to 
phononsiii Such correspondence however does not ap- 
pear to have universal validity. In fact when the ansatz 
(jT5j) is used in the 2D Frohlich- Rashba model, the asymp- 
totic strongly-coupled polaron ground state energy for 
Er — ► 00 does not decreases indefinitely as Eq. (fl"8f but, 
rather, it reaches a minimum finite value.— This is in 
striking contrast with the 3D strongly coupled Frohlich 
polaron in a strong magnetic field, whose ground state 
energy has a squared logarithmic functional form as the 
3D hydrogen atom for H — > 00,— due to the ID confining 
effect of the magnetic field on the electron motion. 
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